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The magnetic excitations of the quasi-2D spin-1/2 frustrated Heisenberg antiferromagnet 
CS2CUCI4 are explored throughout the 2D Brillouin zone using high-resolution time-of- flight inelas- 
tic neutron scattering. Measurements are made both in the magnetically-ordered phase, stabilized 
at low temperatures by the weak inter-layer couplings, as well as in the spin liquid phase above 
the ordering temperature Tjv, when the 2D magnetic layers are decoupled. In the spin liquid phase 
the dynamical correlations are dominated by highly-dispersive excitation continua, a characteristic 
signature of fractionalization of S=l spin waves into pairs of deconfined 5=1/2 spinous and the 
hallmark of a resonating- valence-bond (RVB) state. The boundaries of the excitation continua have 
strong 2D-modulated incommensurate dispersion relations. Upon cooling below Tjv magnetic order 
in an incommensurate spiral forms due to the 2D frustrated couplings. In this phase sharp magnons 
carrying a small part of the total scattering weight are observed at low energies, but the dominant 
continuum scattering which occurs at medium to high energies is essentially unchanged compared 
to the spin liquid phase. Linear spin-wave theory including one- and two-magnon processes can 
describe the sharp magnon excitation, but not the dominant continuum scattering, which instead 
is well described by a parameterized two-spinon cross-section. Those results suggest a cross-over 
in the nature of the excitations from 5=1 spin waves at low energies to deconfined 5=1/2 spinous 
at medium to high energies, which could be understood if CS2CUCI4 was in the close proximity of 
transition between a fractional RVB spin liquid and a magnetically ordered state. A large renor- 
malization factor of the excitation energies (i?=1.63(5)), indicating strong quantum fluctuations in 
the ground state, is obtained using the exchange couplings determined from saturation-field mea- 
surements. We provide an independent consistency check of this quantum renormalization factor 
using measurements of the second moment of the paramagnetic scattering. 

PACS numbers: 75.10.Jm, 75.45.+j, 75.40.Gb, 05.30.Pr 



I. INTRODUCTION 

One of the most remarkable phenomena that can oc- 
cur in strongly correlated systems is the emergence of 
particles with fractional quantum numbers. This occurs 
when fluctuations in the many-body quantum ground 
state are strong enough to provide a screening by which 
fractional components de- confine from local integer con- 
straints. The best known example in magnetism is the 
one-dimensional(lD) S=l/2 Heisenberg antiferromag- 
netic chain (HAFC) where a semi-classical 5'=! spin wave 
spontaneously decays into a pair of deconfined S—1/2 
spinous that separate away indcpendcntly^-^'^. The sig- 
nature of spin fractionalization is a highly-dispersive con- 
tinuum of excited states instead of sharp single-particle 
poles in the dynamical correlations, probed directly by 
inelastic neutron scattering experiments. 

Fractionalization arises from subtle many-body quan- 
tum correlation effects (in the spin-1/2 HAFC chain a 
topological Berry phase term is responsible''^) and higher 
dimensional realizations of such exotic physics have been 
eagerly sought. Amongst the proposed theoretical sce- 
narios is the Resonating- Valence-Bond (RVB) spin liquid 
state introduced by Anderson^ in the context of 2D frus- 
trated quantum magnets. Here spins spontaneously pair 



into singlet bonds which fluctuate between many differ- 
ent configurations; breaking a singlet bond releases two 
5=1/2 spinous that propagate away independently in 
the background of the fluctuating bonds. Central to the 
stability of such an RVB phase is frustration, which en- 
hances quantum fluctuations over the mean-field effects 
that would otherwise favor conventional ordered phases. 

Nearly all experimentally studied phases in 2D quan- 
tum magnets have been found to show conventional con- 
fined phases, which are well characterized. For example, 
in unfrustrated magnets typified by the 5=1/2 Heisen- 
berg antiferromagnet on a square lattice (HSL) mean- 
field effects dominate and quantum fluctuations cause 
only small renormalizations to the semi-classical descrip- 
tion of a Neel ordered ground state with 5=1 transverse 
spin-wave excitations^. Other 2D quantum phases in- 
clude dimerized states recently observed experimentally^ 
where spins are paired into singlet bonds (dimers) ar- 
ranged in a regular (fixed) ordered pattern in the ground 
state; such a phase has a spin gap to a triplet of 5=1 
magnon excitations, i.e. spinous are confined. 

However, spins in frustrated geometries are believed to 
behave in much more unconventional ways, and interest 
in such systems is also motivated by the observation of su- 
perconductivity (possibly mediated by spin fluctuations) 
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in charge-doped triangular lattice materials k-(BEDT- 
TTF)2Xl and Na^jCoOz-yHsO^. In this respect the re- 
cent discovery of a 2D quantum magnet (CS2CUCI4) with 
spins on a triangular lattice which shows fractionaliza- 
tion is of clear importance for elucidating the underly- 
ing magnetismpfi. Indeed, deconfined phases in 2D and 
the conditions required to create them are a major un- 
solved theoretical problem.^iiiii2iiiili^. CsaCuCU is 
a quasi-2D 5=1/2 frustrated Heisenberg antiferromagnet 
on an anisotropic triangular latticaii and neutron scat- 
tering measurements show the dynamical correlations to 
be dominated by a broad continuum of excited statesiS 
as characteristic of deconfined 5=1/2 spinous. Here we 
build on our initial neutron scattering results and pro- 
vide comprehensive measurements throughout the Bril- 
louin zone. Effects of the strong two-dimensionality are 
explicitly observed at all energy scales of the excitations. 
Measurements are made both at temperatures above Ttv 
when the 2D magnetic layers are decoupled as well as be- 
low T/v when mean-field effects from the weak inter-layer 
couplings stabilize 3D magnetic order with an incommen- 
surate spiral structure. 

The paper is organized as follows. The crystal struc- 
ture and magnetism of CS2CUCI4 are described in SecHll 
and the experimental technique used to probe the excita- 
tions is explained in Sec. lIIII The dispersion relation and 
scattering lineshapes measured in the low-temperature 
ordered phase are presented next in Secs. lIVXI and UVBI 
We find sharp magnon peaks carrying a small part of 
the total scattering weight at low energies and highly- 
dispersive continua carrying the majority of the scatter- 
ing weight at moderate to high energies. Results are first 
compared to linear spin-wave theory (reviewed in Ap- 
pendix A) including both one and two-magnon processes 
in Sec. IIVB II This theory is found inadequate to de- 
scribe the dominant continuum scattering, which instead 
is well described by a parameterized two-spinon cross- 
section (Sec. II V B 2|l . Measurements in the spin liquid 
phase above Tjv where the magnetic layers are decoupled 
are shown in Sec. IIV CI The paramagnetic scattering is 
described in Sec IIV Dl where the extracted second mo- 
ment is compared with sum rules. In Sec. IV ll we discuss 
issues in depth with reference to proximity to a spinon 
confinement transition. The paramagnetic scattering is 
discussed in the context of spinon systems at high tem- 
peratures in Sec. IV 2l Finally, the main results and con- 
clusions are summarized in Sec IVII 



II. CRYSTAL STRUCTURE AND MAGNETIC 
PROPERTIES OF CS2CUCL4 

The crystal structure of CS2CUCI4 is 
orthorhombic^^(Fnma) with lattice parameters a— 9.65 
A , &=7.48 A and c=12.26 A at 0.3 K. The structure 
is illustrated in Fig. Qla) and consists of CuCU^" 
tetrahedra arranged in layers {be plane) separated along 
a by Cs"*" ions. The material is an insulator with each 




2 4 6 8 10 12 

Applied field B (T) II a 



FIG. 1: (color online) (a) Crystal structure in CS2CUCI4 
showing the CuCl4~ tetrahedra (pyramids) arranged in lay- 
ers {be plane). The orthorhombic unit cell is indicated by the 
dashed rectangular box. (b) Magnetic exchange paths in a 
{be) layer form a two-dimensional anisotropic triangular lat- 
tice: strong bonds J (heavy lines || b) and smaller frustrating 
zig-zag bonds J' (thin lines), (c) Schematic phase diagram of 
CS2CUCI4 in temperature and magnetic field along a showing 
the region probed by the present experiments (dashed ver- 
tical arrow at B = 0). The magnetic phases are: 3D LRO 
(T < Tjv) with spiral magnetic long-range order, spin-liquid 
{Tm < T < Tmax) characterized by strong intra-layer an- 
tiferromagnetic correlations (at r,„ai=2.65 K the magnetic 
susceptibility has a maximum), paramagnetic (T 3> (J, J')) 
and ferromagnetic {B > Be) where spins are ferromagneti- 
cally aligned by the applied field. The solid line is a phase 
transition boundary and dashed lines show cross-overs. 



Cu^"*" ion carrying a spin of 1/2. Crystal field effects 
quench the orbital angular momentum resulting in 
near-isotropic Heisenberg spins on each Cu^^ ion. There 
are four such ions in each orthorhombic unit cell, two 
located on each CuCl layer as illustrated in Fig. ^a) . 

Our previous measurements^^ showed that CS2CUCI4 
is a quasi-2D low-exchange quantum magnet. This quasi- 
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2D character is a result of the layered crystal structure, 
which restricts the main superexchange routes to neigh- 
boring spin sites in the (be) plane. Since this superex- 
change route is mediated by two nonmagnetic Cl~ ions 
the exchange energies are low ~l-4 K. In each layer 
the exchange paths form a triangular lattice with non- 
equivalent couplings indicated in Fig. ^b): exchange J 
along the &-axis and J' along the zig-zag bonds. The 
main 2D Hamiltonian describing CS2CUCI4 is thus 

W= ^ JS,-5,. + ^ -Sj, (1) 

with each interacting spin-pair counted once. This 
Hamiltonian interpolates between non-interacting chains 
(J' — 0), the fully frustrated triangular lattice (J' — J), 
and unfrustrated square lattice (J = 0). Since the main 
spin couplings define a 2D isosceles triangular lattice as 
shown in Fig. ^b), measurements of the excitations will 
be discussed in terms of the 2D Brillouin zone of this tri- 
angular lattice, even though the full crystal symmetry is 
orthorhombic. 



A. Hamiltonian 

The full spin Hamiltonian of CS2CUCI4 and its pa- 
rameters have been previously determined^^ from mea- 
surements in high applied magnetic fields that overcome 
the antiferromagnetic couplings and stabilize the fully- 
aligned ferromagnetic ground state | ttt ■ • ■) in Fig. 
^c). In this unique phase quantum fluctuations are en- 
tirely suppressed by the large field, the excitations are 
magnons and their dispersion relations image directly 
the Fourier transform of the bare (i.e. unrenormalized) 
exchange couplings. Using this method the determined 
interactions are: J=0. 374(5) meV along the &-axis and 
J'/J=0.34(3) along the zig-zag bonds, the couplings in 
Fig. ^b), and small Dzyaloshinskii-Moriya terms along 
the zig-zag bonds Da/ J— 0-053{5) and inter-layer cou- 
plings J"/J=0.045(5) (for further details see Ref J ig). 
Note that sum rules to be discussed in Sec. IIVDI also 
provide an independent verification of the strength of the 
couplings. 

The finite interlayer couplings (J") stabilize long-range 
magnetic order— below Tjv=0.62(l) K. The order is in- 
commensurate due to the frustrated couplings in the (be) 
plane with the ordering wavevector Q — (0.5 -I- eo)b*, 
£0=0.030(2). Ordered spins rotate in a spiral nearly con- 
tained in the (be) plane due to the small DM couplings, 
which create an effective easy-plane anisotropy. The pro- 
jection of the spiral order onto the (be) plane is schemat- 
ically illustrated in Fig. lTUT b'). 

III. EXPERIMENTAL DETAILS 

The samples used were high-quality, large single crys- 
tals of CS2CUCI4 grown from solution^- and cut into 



plates of typical size 2 x 2 x 0.6 cm^ optimized for neutron 
absorption. Temperature control was provided by a dilu- 
tion refrigerator insert with a 0.04 K base temperature. 
Measurements were made in the different regions of the 
magnetic phase diagram indicated in Fig.^c) (zero-field 
line): the 3D spiral ordered phase (below 0.1 K), the spin 
liquid phase above Tjv (at 0.75 and 0.9 K) and also the 
paramagnetic phase (at 12.8 and 15 K). 

The magnetic excitations were measured using the 
indirect-geometry time-of-flight spectrometer IRIS^^ at 
the ISIS spallation neutron source in the UK, which 
combines high energy resolution with a large detector 
solid angle that allows simultaneous coverage of extended 
(wavevector, energy) regions. The detector array collects 
data on a near-cylindrical surface in (k^ ,ky,E) space 

(kx - \/[Ef -\- E)/a^ +ky — Ef/a, where k^ (along the 

incident beam direction) and ky are the components of 
the wavevector transfer in the horizontal scattering plane, 
E is the energy transfer, Ef is the fixed final energy and 
a = 5.072 meV A^. Dispersion points are obtained when 
this measurement surface intersects the dispersion sur- 
face of the magnetic excitations (a typical intersection is 
shown by the bold curve (2) in Fig. |3a)). By changing 
the orientation of the crystal the measurement surface 
can be positioned to probe the excitations near particu- 
lar directions in reciprocal space. Several non-equivalent 
directions (curves) in the Brillouin zone were probed and 
are indicated by the dotted lines numbered (1) to (4) in 
Fig.EIb). Typical counting times for one crystal orienta- 
tion and temperature ranged between 30 to 54 hours at 
an average ISIS proton current of 170 /iA. Capital letters 
A-K indicate location of energy scans with parameters 
listed in Table. Throughout this paper the wavevec- 
tor transfer k = (h, fc, I) is expressed in reciprocal lattice 
units (rlu)22 of (27r/a, 27r/6, 2it/c). 

An energy resolution Ai?=0. 019(1) meV (fuU-width- 
at-half-maximum)i^ on the elastic line was achieved us- 
ing a fixed final energy £'/ = 1.846 meV selected by py- 
rolytic graphite (002) analyzers, cooled to liquid Helium 
temperatures to reduce the thermal diffuse background. 
The detector array had 51-elements covering the range 
of scattering angles 25.75° < 29 < 158.0°. The raw 
neutron counts vs. time-of-flight were converted into 
scattering cross-section intensities S{k,Lo) [typical data 
is shown in Fig. 0] , which were further corrected for neu- 
tron absorption effects using a numerical calculation for 
a plate-shaped crystal calibrated against the measured 
incoherent quasi-elastic scattering from the sample. The 
non-magnetic background was estimated from measured 
intensities near places in the Brillouin zone such as the 
ferromagnetic zone center where the magnetic scatter- 
ing contribution is negligible; the accuracy of this back- 
ground subtraction procedure is illustrated in Fig. [SfJ 
(solid points). 
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TABLE I: Scan parameters: wavevector fc — [h, k, I) and energy [E) change along scan direction, location of main peak Uk and 
polarization factors {px,Pz) at this point (see Eq. ljA12fl '). 



Scan 


/i(rlu) 




fc(rlu) 




/(rlu) 


£(meV) 


cjfe(meV) 


Px 


Pz 


A 


1.31 + 0.29-E - 


0.025£;^ 


-0.389 + 0.189S 


- 0.016^;^ 





E 


0.69(2) 


1.95 


0.05 


B 


1.69 + 0.29E'- 


0.025£;^ 


-0.30 + 0.189E' 


- 0.015^^ 





E 


0.62(2) 


1.98 


0.02 


C 


0.84 




0.21 + 0.297B - 


- 0.026£;^ 





E 


0.56(2) 


1.75 


0.25 


D 


0.69 




2.11 + 0.29S - 


0.025^2 





E 


0.68(2) 


1.05 


0.95 


E 







-0.33 + 0.19S- 


- 0.015^^ 


0.78 + 0.37B - 0.03£;^ 


E 


0.56(1) 


1 


1 


F 







-0.39 + 0.19^; 


- 0.02£^ 


1.66 + 0.37S- 0.035£;^ 


E 


0.67(2) 


1 


1 


G 







0.5 




1.53 - 0.32B - Q.IE^ 


E 


0.107(10) 


1 


1 


H 







0.28 + 0.29S - 


0.025£^ 


1.205 


E 


0.34(2) 


1 


1 


I 







k 




-0.23 + 3.63fc - 1.64fc2 


0.85(5) 




1 


1 


J 







0.47 




1.0 - 0.45B 


E 


0.10(1) 


1 


1 


K 







0.29 + 0.29S- 


- 0.03£;^ 


0.77-0.14£ + 0.013£;^ 


E 


0.32(2) 


1 


1 



TABLE IL Dispersion parameters. 



Zone-boundary energy (meV) along [OfcO] 

Zone-boundary energy (meV) along [OOZ] 

measured at fe = (0,0.5, 1.5) 
Dispersion (meV) along / for k ~ —0.25 
Dispersion (meV) along h for k ~ —0.25 
Incommensuration eo(rlu) 



0.67(1) and 0.55(1) 

alternating 

0.107(10) 

0.12(1) 
< 0.02 
0.030(2) 



IV. RESULTS AND ANALYSIS 

A. Dispersion relations in the spiral ordered phase 
below Tjv 

Measurements of the magnetic excitations at the low- 
est temperature (T < 0.1 K) in the spiral ordered phase 
showed strong dispersion along both directions in the 2D 
triangular layers ([OfcO] and [00/]) and negligible disper- 
sion along the inter-layer direction [/lOO] . Table |n] lists 
key parameters and Fig. [3 shows a surface plot of the 
dispersion relation in the 2D triangular plane. Data ob- 
tained from cuts through this dispersion surface along 
several non-equivalent directions in the 2D plane [curves 
numbered (1) to (4)] is plotted in Fig. |31 

Fig- Ell) gives the dispersion along fc for /=0. Here 
data points also have a finite component along h, but 
no dispersion could be detected along this direction for 
near-constant fc, as expected for weak couplings J" be- 
tween layers. The dispersion along fc is incommensu- 
rate, with the minimum energy at the magnetic Bragg 
peak positions indexed by the ordering wavevector Q — 
(0.5 -I- eo)b*, £0=0.030(2), and has an asymmetric shape 
with non-equivalent zone boundary energies of 0.67(1) 
and 0.55(1) meV on the two sides of the minimum. The 



incommensuration and asymmetry of the dispersion are 
direct consequences of the frustrated couplings (a sym- 
metric dispersion, commensurate at the antiferromag- 
netic point fc=0.5 would be expected in the absence of 
frustration J'=0 in Fig.IJb)). The effects of the 2D cou- 
plings are most clearly seen in the dispersion along I for 
near-constant fc ^ —0.25 plotted in Fig. |3j2). The ob- 
served dispersion is 0.12(1) meV, practically equal to the 
strength of the couphngs along I, J'=0. 128(5) mcV as 
directly measured in Ref. ^3- The significant dispersion 
observed along both [OfcO] and [00/] directions indicates 
a strong 2D character for the magnetic excitations. 

Further measurements that illustrate the /-dependence 
of the dispersion are shown in Fig.|3l3), with data near 
/ = 1. Compared to /=0 data in (1) the left-right asym- 
metry of the dispersion along fc has now been reversed 
and the minimum energy displaced to the new Bragg 
peak position in this zone fc = 0.5 — eq, as expected by 
periodicity for a triangular lattice with a 2D Brillouin 
zone as shown in Fig. [3 Data near the zone-boundary 
point along /, /=1.5, are plotted in Fig.|3f4). Here the dis- 
persion along fc appears nearly symmetric around fc=0.5 
with a gap at W(o^o. 5, i.5)=0. 107(10) meV. 

The observed modulations in the dispersion relation as 
a function of fc and / can be well described by the prin- 
cipal spin-wave dispersion of the 2D frustrated Hamilto- 
nian in Eq. (for details see Appendix A) 



{Jk - Jq) [{Jk-Q + Jk+q) /2 - Jq] (2) 



where the Fourier transform of the exchange couplings 
is Jfc — Jcos27rfc -|- 2J'cos7rfccos7r/, k = {h,k,l)^ and 
the ordering wavevector Q = (0.5-|-eo)b* as above. Solid 
lines in Fig. El show the intersection between the scan 
directions and Eq. ||2J) (using the actual fc, / and E at 
each point in the scan). The effective exchange param- 
eters obtained from the above global fit are J=0.61(l) 
meV and J'=0. 107(10) meV, in agreement with earlier 
estimatesia (J=0.374(5) meV and J'=0. 128(5) meV are 
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k ilTilb) 

FIG. 2: (color online) (a) Surface plot of the dispersion re- 
lation in Eq. Q as a function k and I with constant-energy 
contours shown projected onto the basal plane (thin curves), 
(b) Reciprocal space diagram showing where measurements 
were made. Numbers (1) to (4) refer to directions (curves) in 
reciprocal space where dispersion points (solid dots) were ex- 
tracted (curves (1) and (2) are also indicated by bold curves 
on the surface plot in (a)); the extracted dispersions are plot- 
ted in the numbered panels in Fig. |3 Capital letters A-J 
indicate location of energy scans. Thick lines show the near- 
hexagonal Brillouin zones of the 2D triangular lattice in Fig. 
^b), filled circles are zone centers (t) and stars at r ± Q, 
Q = (0.5-|-eo)b* mark incommensurate magnetic Bragg peaks 
from the spiral order at T < Tjv. 



the hare exchange energies measured using the high-field 
technique ^2, where quantum fluctuations are quenched 
out). The renormalization of the excitation energy com- 
pared to the classical spin-wave result J/J=1.63(5), is 
very large and is similar to the exact result 7r/2 for the 
ID 3=1/2 HAFC (see e.g. Ref. 0), where it originates 
from strong spin-singlet correlations in the ground state. 
In contrast, the spin- wave velocity (energy) renormaliza- 
tion in the unfrustrated S=l/2 HSL is only 1.18, see Ref. 
0. Such a large renormalization indicates strong quan- 
tum fluctuations in the ground state of CS2CUCI4. 

The Z-dependence of the low-energy dispersion is il- 
lustrated in Fig. 1215), which plots the oscillation in the 



minimum energy position vs. k and /. Data points were 
extracted form flts to the low- to medium-energy part 
(0.1< E <0.4 mcV) of the dispersion along k for various 
Vs. The observed incommensurate oscillation is well ex- 
plained by Eq. (0), for which at a given I the minimum 

energy occurs at fc = 0.5+^ sin^^ ( ~ costtI) (solid line). 



B. Excitation lineshapes in the spiral ordered 
phase below Tjv 

Lineshapes of scattering continua are important as 
they give information about the underlying quasi- 
particles and their interactions. In CS2CUCI4 the dynam- 
ical correlation are dominated by extended continua as 
indicated by the dotted areas in Fig.O A representative 
scan near the antiferromagnetic zone boundary along k is 
shown in Fig. ^bottom). The scattering is highly asym- 
metric with a significant high-energy tail. Instrumental 
resolution effects are minimal as the overall extent of the 
inelastic scattering is an order of magnitude larger than 
the instrumental resolution (grey shaded peak). The 
non-magnetic background (dashed line) is modelled by a 
constant-plus-exponential function. The magnetic peak 
disappears at 15 K (solid circles) and is replaced by a 
broad, overdamped, paramagnetic signal. Instead of con- 
taining sharp peaks, characteristic of conventional, 5=1 
magnon excitations, the observed low-temperature line- 
shape is dominated by a broad continuum scattering. To 
quantify the extent of the magnetic scattering we de- 
fine the "upper boundary" u!^ as the energy above which 
the observed intensities could not be distinguished from 
the non-magnetic background. The upper boundary dis- 
persion is shown in Fig. |21 (large open circles and upper 
dashed line) and is included to indicate the region below 
which most of the magnetic scattering is located. 

Representative scans throughout the Brillouin zone are 
shown in Fig. where in this case the data is properly 
normalized and corrected for absorption effects, and with 
the non-magnetic background subtracted. Sharp thresh- 
olds observed at the lower boundaries of scattering con- 
tinua, most clearly seen in scans A and F, are indica- 
tive of highly-coherent excited states, as opposed to the 
broad, featureless scattering in systems with spin freezing 
or random disorder. 

We first analyze the scattering lineshapes in terms of 
linear spin- wave theory (LSWT), which provides a good 
description of the excitations in unfrustrated 2D square- 
lattice Heisenberg antiferromagneta^. 



1. Lineshape analysis using linear spin-wave theory 

The spin-wave theory for a frustrated Hamiltonian 
with a spiral-ordered ground state (reviewed in Appendix 
A) predicts a principal spin-wave mode cok polarized per- 
pendicular to the plane of spin rotation (out-of-plane) 
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FIG. 3: Dispersion relation of the magnetic excitations at T <0.1 K. Numbered panels (l)-(4) correspond to directions in the 
2D plane indicated in Fig.|51;b) [(1) is from Ref. ,10j. Filled symbols are the main peak in the lineshape and the solid line is a 
fit to the spin- wave dispersion relation in Eq. ((2J; the dotted area indicates the extent of the magnetic scattering, large open 
circles mark the experimentally-estimated continuum upper boundary (the upper thick dashed line is a guide to the eye) and 
open squares show the lower boundary of the scattering (where it could be resolved from the main peak as shown in Fig. 0J ; 
dashed and dash-dotted lines are the secondary spin- wave dispersions and , respectively; light shaded regions labelled 
with capital letters A-K (the line thickness represents the wavevector averaging) indicate scan directions. Panel (5) shows 
the incommensurate modulation in the dispersion relation vs. k and I compared with predictions of Eq. @ (solid line) and 
uncoupled chains ( J' = dashed line). 



and two secondary modes = Wfc+Q and uj'j^ = ujk-Q, 
polarized in-plane. The three dispersions are graphed in 
Fig. 13 and the corresponding excitations can be identi- 
fied with quantized modes where the total spin compo- 
nent along the z-axis normal to the spiral plane changes 
by AS''=0 (wfc), +1 (w^) and -1 (u;^) {+z defines the 
sense of rotation in the spiral) . 

(a) Sharp magnon excitation at low energies 

Scans that probe the low-energy excitations observe a 



sharp peak below the onset of the higher-energy contin- 
uum scattering. Fig. [SfJ shows measurements near the 
zone-boundary point along I, fe=(0,0.5,1.5) [scan G in 
Fig. El3)]. A sharp peak occurs at Wfe^O. 107(10) meV 
followed by strong continuum scattering with an onset 
at a slightly higher energy around uj* ~0.2 meV and ex- 
tending up to 1 meV. The sharp low-energy peak can 
be well described by the spin-wave cross-section in Eq. 
HA10|I (dashed line), which predicts a dominant (out-of- 
plane polarized) magnon at the lowest energy followed by 
very much weaker (~ 1/6 intensity) secondary modes. 
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0.0 0.2 0.4 0.6 0.8 1.0 1.2 
Energy (meV) 

FIG. 4: (Bottom) Intensity observed in scan D in the spiral 
ordered phase at 0.1 K (open symbols) and in the param- 
agnetic phase at 15 K (filled circles). Data points are raw 
neutron counts and solid lines are a guide to the eye. The 
dashed line is the estimated non-magnetic background and 
the grey shaded peak indicates the instrumental energy res- 
olution. (Top) Intersection of scan D with the dispersion re- 
lations. Vertical arrows between the top and bottom panels 
show the onset of the magnetic scattering (open square), the 
main peak in the lineshape (solid circle), and the upper con- 
tinuum boundary (open circle), respectively. 



(b) Continuum scattering vs. two-magnon processes 

The largest intensity in scan G is however contained in 
the continuum scattering (about 2/3 of the total energy- 
integrated intensity). In contrast, LSWT predicts a much 
smaller continuum intensity arising from two-magnon 
scattering (about 15%). In fact the shaded area in Fig. 
I^JGI shows the two-magnon intensity scaled upwards by a 
factor of 9 to illustrate it (the calculation is described in 
Appendix A and is based on Eq. ljAll|l and sum rules for 
the total scattering). The two-magnon lineshape would 
predict the onset of the continuum scattering immedi- 
ately above the one-magnon energy (jjfc=0. 107(10) meV, 
whereas a distinct separation is observed in the data 
with the continuum scattering starting around u* ~0.2 
meV. The predicted two-magnon intensity drops off very 
fast with increasing energy (especially at high energies 
E >0.6 meV~ J), whereas the observed scattering falls 
off significantly slower [see inset in (Spl] . 

The measured continuum scattering as a function of 
wavevector at constant energy (scan I in Fig. 12^3)) is 
compared against the two-magnon cross-section in Fig. 



O. This energy is beyond the one-magnon cutoff and 
only two-magnon processes contribute, and their inten- 
sity is shown by the shaded area(the same x9 scaling 
factor as for scan G was used). Apart from largely under- 
estimating the overall strength of the observed scatter- 
ing intensity at this energy, the two-magnon functional 
form is also very different from the data, predicting a 
slightly wider extent in wavevector and intensity modula- 
tions (peak at k ~0.75 originating from an increased den- 
sity of states for two-magnon scattering) that are not ob- 
served in the data. Similar disagreement occurs between 
the observed continuum scattering and a two-magnon 
lineshape at other wavevectors throughout the Brillouin 
zone. From this analysis we conclude that the observed 
functional form (energy- and wavevector-dependence) of 
the observed continuum scattering is not reproduced by 
a two-magnon cross-section. 

The sharp peak observed below the scattering contin- 
uum in scans such as in Fig. [SfJ could only be clearly re- 
solved at low energies below 0.2 meV. With increasing 
Wfe this sharp peak appears to gradually merge into the 
higher-energy continuum scattering and a typical line- 
shape is shown in Fig. |Sfl [scan H in Fig. El3)]. The 
intensity observed here is dominated by continuum scat- 
tering with a power-law lineshape above a lower bound- 
ary Wfc=0.34(2) meV (solid line, Eq. Q) and no sepa- 
rate sharp mode could be resolved below the continuum 
lower boundary. A spin-wave cross-section Eq. (jA10|) 
(dashed line) scaled to the intensity of the one-magnon 
peak in scan G, significantly underestimates the inten- 
sity observed near the lower boundary of the continuum 
scattering. At higher energies the spin- wave cross-section 
is also inadequate to describe the low-energy part of the 
scattering lineshapes as shown by dashed lines in Fig. 
|3K,F. 

(c) Polarization of the continuum scattering 

In principle the polarization of the modes can be ex- 
tracted from neutron scattering data because of the di- 
rectional dependence of the intensities. We now consider 
this property to find the polarization of the scattering 
with respect to the plane of the spiral order. We com- 
pare the intensities measured at two equivalent places in 
the 2D Brillouin zone [A and F in Fig. |2Ib)], but where 
the relative weights of the in-plane and out-of-plane po- 
larizations observed by neutron scattering are very dif- 
ferent. In Fig. |S1 scan A observes mainly in-plane polar- 
ized scattering, whereas scan F is equally sensitive to in- 
plane- and out-of-plane polarizations, see Tabled Con- 
trary to the polarized spin-wave cross-section (dashed 
lines, Eq. (|A10|) ). which predicts significant changes in 
the relative intensity of the different components of the 
lineshape due to the changes in the polarization factors, 
the two scans observe very similar lineshapes and the 
scattering is instead consistent with an isotropic polar- 
ization. The high-energy continuum scattering is essen- 
tially unchanged upon heating above T/y (see Fig. ^JGI-I) , 
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FIG. 5: Magnetic inelastic scattering measured in the ordered phase below 0.1 K in scans A- J. Solid lines are fits to a modified 
two-spinon cross-section, Eq. (|1J for A-D and Eq. Q for G-J. Vertical arrows indicate the estimated continuum upper boundary 
io^ . Dashed lines show predicted lineshape for polarized cycloidal spin waves Eq. HAIOII and the dark shaded region (D,G,I) 
indicates the predicted two-magnon scattering continuum calculated using Eq. lIXTil and sum rules (scaled up by a factor of 9 
in G and I) . All calculations include the convolution with the spectrometer resolution function and the isotropic magnetic form 
factor for Cu^^ ions^^ . Solid points in G show the accuracy of the background subtraction for wavevectors where no magnetic 
scattering is expected, squares for fc=0.70(5), E < 0.4 meV, and filled circles for fc=1.00(5), E >0.4 meV in Fig.0l4). 



also consistent with a mainly isotropic polarization of the 
continuum scattering, independent of the orientation of 
the spiral ordering plane. 

(d) Sharp mode near saddle-points 

Sharp peaks were observed at high energies near special 
wavevectors where the 2D dispersion uj^ is at a "saddle" 
point as indicated in Fig. HJa) (such points are located 
near fe=(0,±0.25,/) with I an odd integer and equiva- 
lent positionsSSi). A typical scan is shown in Fig. 
The lineshape here is dominated by a sharp, resolution- 
limited peak at ti^(o,-o.23,o.98)=0-56(l) meV, followed by 
only very weak continuum scattering at higher energies. 
This sharp mode is much more intense (by a factor of 
4.4(8)) than predicted by extrapolating to those ener- 
gies the intensity of the low-energy magnon observed 
at 0.107(10) meV in scan G assuming a spin-wave de- 
pendence, Eq. (jA7|l . suggesting an anomalous intensity 
vs. wavevector dependence or different origin (perhaps 



a multi-particle bound state) for the saddle-point mode. 
Also, spin- wave theory predicts an additional mode near 
a;^=0.7 meV (dashed line), which is however not ob- 
served in the data. The saddle-point mode is consistent 
with an out-of-plane polarization by comparing its in- 
tensity to that of the isotropic continuum scattering ob- 
served near fc=0.5. 

From this analysis we conclude that LSWT is inade- 
quate to describe the results, in particular the dominant 
continuum scattering, and next we consider another ap- 
proach based on experience in ID systems where LSWT 
breaks down and highly-dispersive continua are observed. 

2. Lineshape analysis using a two-spinon cross-section 

Neutrons scatter by changing the total spin by ASt — 
0,±1 and the observation of dominant scattering con- 
tinua indicates that the elementary quasi-particles cre- 
ated in such processes are not conventional S=l spin 



waves. A natural explanation is that quasi-particles have 
fractional spin quantum number (spin is a good quantum 
number since the main (J, J') Hamiltonian has isotropic 
exchanges). By analogy with the ID 5*= 1/2 HAFC chain 
[J' — 0] we identify those quasiparticles with spin- 1/2 
spinous, created in pairs in a neutron scattering process. 
The dominant neutron scattering cross-section is then a 
two-spinon continuum, isotropic in spin space and ex- 
tending between dispersive lower and upper boundaries 
obtained by convolving two spinon dispersions. 

We compare the observed continuum scattering with a 
generic power-law lineshape 



S{k,u)^I, 



e{uj- 



CO) 



,2ll-,,/2 



(3) 



where ujk and uj^ are the experimentally- determined 2D- 
dispersive lower and upper boundaries of the continuum 
scattering, Ik is a wavevector-dependent intensity fac- 
tor. This form is generalized form the two-spinon cross- 
section in the ID HAFC chain (?7=1) and has also been 
proposed in a theoretical description of 2D frustrated 
quantum magnets with deconfined spinous^. 

Scans that probe the continuum scattering over a wide 
energy range are shown in Fig. [3~J,J. The intensities ob- 
served above a cross-over energy scale cu* ~0.2 meV can 
be well described by a power-law lineshape and fits (solid 
lines) to Eq. ^ (assuming lower boundary fixed by the 
dispersion Wfc) give an exponent 77=0.74 ± 0.14 (the inset 
in Fig. shows a log plot of the continuum scattering 
in scans G and J and the solid lines are power-law fits). 
The observed wavevector dependence of the continuum 
scattering at constant energy is also consistent with the 
functional form in Eq. (|3J) as shown by the solid line in 
Fig.E. 

At places close to the zone boundary such as Fig. [Sf 
the observed lineshape can be described by a superpo- 
sition of two continua, each with a power-law lineshape 
as in Eq. Q and onsets at separate lower boundaries. 
In the figure the solid line shows a superposition of two 
continua with onsets at a;^=0.53(2) meV and ti'fe=0.67(2) 
meV. The scattering at 1=0 shown in Fig. |S|^-D can be 
described by a superposition of three continua with lower 
boundaries at the dispersion relations uik, and w^, 
isotropic in spin space and with equal relative weights 
(solid lines): 



6{UJ - UJk) 



.2)l-')/2 (^2 _ (^^ )2)l-./2 



(^2_(^+)2)l-,/2 



(4) 



where the exponent rj has values in the range 0.7 to 1. 

To summarize, throughout most of the Brillouin zone 
and over the largest energy scales probed the dynami- 
cal correlations can be well-described by continuum line- 
shapes, characteristic of a two-spinon cross-section. 
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FIG. 6: Magnetic inelastic scattering measured at 0.75 K in 
the spin liquid phase above Tjv (open symbols) and at 0.05 K 
below Tjv (solid symbols). Top axis shows wavevector change 
along scan direction. Solid lines in E-F are fits to a damped- 
harmonic oscillator lineshape Eq. @ and in G-I are guides to 
the eye (dashed line in G shows a fit to the one-magnon cross- 
section Eq. lAlOt . Data very close to the elastic line E—0 was 
omitted in G. 



C. Dispersion relation and excitation lineshapes in 
the 2D spin-liquid phase above Tjv 

The excitations are thermally broadened which results 
in a smooth low-energy tail for the scattering lineshapes 
as shown in Fig. |Sp by data at 0.75 K in the spin liq- 
uid phase above Tjy. The sharp magnon mode observed 
at base temperature below the continuum lower bound- 
ary in Fig.|SlGI(solid points) at tjfc=0. 107(10) meV damps 
out upon heating (open symbols) and merges with the 
low-energy tail of the continuum scattering. At high en- 
ergies, however, the continuum scattering is practically 
unchanged by small temperatures, see Fig. 

The dispersion relation above T/v maintains a strong 
2D character. The observed incommensuration in the 
low energy part of the dispersion relation is plotted as 
a function of k and I in Fig. |3f5)(open symbols) and 
is essentially unchanged compared to base temperatures 
(solid symbols). The excitations disperse strongly along 
both directions in the 2D plane and the observed disper- 
sion along I for near-constant k —0.25 is plotted in 
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FIG. 7: Dispersion along I for near-constant k ~-0.25 mea- 
sured at 0.75 K in the spin liquid phase above Tjv. The solid 
line is Eq. (|5J with the same parameters as below Tm- Inten- 
sities observed along scans E-F are shown in Fig.|S] 



Fig. [T^solid symbols) and is practically the same as be- 
low (solid line) . The dispersion points were extracted 
by fitting the observed lineshapes to a damped-harmonic 
oscillator form (difference of two Lorentzians centered at 

LO = ±UJk) 



S{k, 



(w2 



[n{uj) + l], (5) 



where ujk is the dispersion relation, is the damping 
rate and n(uj) — l/[exp{Lo/kBT) — 1] is the Bose popu- 
lation factor. This form gave a good description of the 
data as shown by typical fits (solid lines) in Fig.Efl-F. 



D. Paramagnetic scattering 

Upon further increasing the temperature the lower 
boundary collapses and in the paramagnetic phase in- 
elastic intensity is observed in the form of a broad con- 
tinuum scattering extending up to a dispersive upper 
boundary. Fig. IS^top left) indicates the extent of the 
scattering (dotted area) observed at 12.8 K{—2.9 J) as 
a function of wavevector and energy. The intensity ob- 
served in scan K near the antiferromagnetic zone center 
along [OfcO](fc '--^O.S) where the scattering has the largest 
extent in energy is shown in Fig. |SJbottom). The para- 
magnetic lineshape (open circles) is broad with scatter- 
ing extending up to nearly the same upper boundary as 
the continuum scattering observed at low temperatures 
(filled circles). 

We quantify the energy width of the scattering 
in terms of the normalized second moment (ujI) — 
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FIG. 8: (Top left) Extent of paramagnetic scattering observed 
at 12.8 K (dotted area). The same wavevector-energy space 
was probed as for the 0.05 K data in Fig. |32). Large open 
circles are the upper boundary of the scattering with the 
thick (upper) dashed line a guide to the eye. Solid points 
are the square root of the second moment of the scattering 
{luD compared to predictions based on sum rules (solid 
line Eq. ®); thin dashed line shows the same calculation for 
T = CO. (Top right) Intensity observed along scan K. Top axis 
shows wavevector change along scan direction. Thick dashed 
line is a fit to a plateau-like function and dash-dotted line 
shows a fit to a Gaussian with the same second moment. Ver- 
tical arrow indicates the estimated upper boundary of the 
scattering. (Bottom) Intensity along scan K in the paramag- 
netic phase (open circles 12.8 K) compared to low tempera- 
tures (solid circles 0.05 K). Dotted region highlights the area 
of the paramagnetic scattering. Solid line is a fit of the low- 
temperature data to a power-law lineshape Eq. 



J du) Lo^ S{k^uj) / J duj S{k,uj), which can be directly re- 
lated to the exchange couplings in the spin Hamiltonian 
using sum rules. Following Refs . l27l and l28l we obtain for 
a Heisenberg Hamiltonian to first order in f3 — l/ksT 

i^l) = ls{S+l)J2Js{^~cos{k-S)]x 



(6) 



where the sum is over all exchange couplings Jg, S is a 
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vector connecting sites and Jfc = 5 Js exp(ifc-5) is the 
Fourier transform of the exchange couphngs. For the ID 
5=1/2 HAFC model (J'=0 in Eq. ((TJ) the above result 
predicts that the paramagnetic scattering ai T = 00 has 
a sinusoidal width \/Juj^ — \/2>/| sin7rfc| that is zero 
at the ferromagnetic zone center and increases up to a 
maximum of \/2J at the antiferromagnetic zone center. 

The second moment was extracted from scans such as 
in Fig.lHl^top right) where the data on the negative energy 
side was constructed from the data on the positive energy 
side according to the principle of detailed balance. This 
procedure gave estimates of (w^) averaged over a small 
range of fc-values and the results are plotted in Fig.lSJtop 
left) (solid symbols). The observed second moment near 
the antiferromagnetic zone center and its reduction upon 
approaching the ferromagnetic zone center are well re- 
produced by Eq. © (solid line) for the 2D Hamiltonian 
in Eq. ^ with the bare exchange couplings J =0.374 
meV and J'=0.128 meV as determined for CS2CUCI4 us- 
ing saturation-field measurementsii. Although this high- 
temperature method of probing the Hamiltonian does not 
give the exchange couplings J and J' directly, it does 
however provide a measure of the overall energy scale 
of the couplings without involving magnetic fields and is 
thus an independent consistency check of the bare ex- 
change couplings and quantum renormalization factors. 
The second-moment calculation predicts that no signifi- 
cant changes are expected upon heating to T = 00 (thin 
dashed line in Fig. |SJtop left)) showing that the mea- 
surement temperature (2.9J) was sufficiently deep into 
the paramagnetic phase. 



V. DISCUSSION 

1. Magnetic excitations and cross-over phase diagram 

Our experiments on CS2CUCI4 described above observe 
that in the spin liquid phase above Tjv the dynamical 
correlations are dominated by highly-dispersive scatter- 
ing continua, the hallmark of deconfined 3=1/2 spinous 
of a fractionalized phase. Upon cooling in the spiral 
ordered phase below T^r an 5=1 magnon with a small 
scattering weight emerges at low energies as expected 
for a phase with a continuous broken symmetry, how- 
ever the dominant scattering still occurs in the form of 
a continuum of excitations at medium to higher ener- 
gies and is best described in terms of pairs of deconfined 
spinous. This could be understood if CS2CUCI4 was in 
the close proximity of a fractionalized phase with decon- 
fined spinous, but only weakly perturbed by small terms 
in the Hamiltonian. Those weak terms stabilize long- 
range magnetic order and create an effective short-range 
attractive potential between spinous that is sufficient to 
create a two-spinon bound state (5=1 magnon) at low 
energies, but deconfined spinous still occur at energies 
higher than the small scale of the attractive potential. In 
the spin liquid phase above small thermal fluctuations 



overcome this small perturbation in the Hamiltonian and 
only deconfined spinous would be expected, in agreement 
with experiments which observed two-spinon continua at 
all energy scales probed. Proximity of CS2CUCI4 to a 
spinon-binding transition is supported by the observa- 
tion of sharp divergencies at the lower boundary of two- 
spinon continua (most clearly seen in scans A and F), 
which indicate a resonant enhancement of the scattering 
as precursor to the formation of a sharp bound state, as 
noted in Ref. Moreover, evidence that the ordered 
state is stabilized by only small terms in the Hamilto- 
nian comes form measurements in in-plane fields which 
observe suppression of the spiral order and transition to 
a disordered phase at relatively small applied fieldsiS. 

An important question is the nature of the fractional- 
ized phase that dominates the physics in CS2CUCI4 and 
leads to strong two-spinon continua in the dynamical cor- 
relations. Various approaches have been actively pur- 
sued theoretically. A quasi- ID approach was proposed^, 
starting from ID HAFC chains along the strongest- 
exchange direction in the 2D triangular planes and treat- 
ing perturbatively all other terms in the Hamiltonian, 
including the frustrated coupling J'; in this scenario the 
relevant picture is weakly-coupled chains'^^ where spinous 
are ID objects identified with quantum solitons. How- 
ever, in CS2CUCI4 the coupling between "chains" is rel- 
atively large J' = 0.34(3) J and moreover, the measured 
dispersion relations show strong (of order J') 2D modula- 
tions at all energy scales, indicating that the excitations 
are strongly affected by the two-dimensionality of the 
couplings. Explicit 2D approaches to the physics, such as 
the resonating- valence-bond (RVB) picture^, have been 
considered by other authors^^-^'^. The RVB state can be 
phenomenologically described as a spin liquid where spins 
are spontaneously paired into singlet bonds that fluctu- 
ate (resonate) between many different conflgurations to 
gain quantum kinetic energy. In this phase spinous are 
2D objects with topological charactesi^ and can be physi- 
cally described as the two spin- 1/2 ends of a broken bond 
that separate away through bond rearrangement. Ex- 
plicit calculations within the two approaches (quasi-lD or 
explicit 2D) are required to compare with the measured 
dispersion relations and scattering lineshapes to deter- 
mine which one best captures the physics in CS2CUCI4. 

If the fractionalized phase relevant to CS2CUCI4 is a 
2D RVB phase then a possible phase diagram scenario 
showing the transition to incommensurate 3D spiral or- 
der driven by some small parameter in the Hamiltonian 
(such as the interlayer couplings J") is illustrated in Fig. 
|51 This phase diagram is inspired from the proposal^^' — 
that if a 2D frustrated quantum magnet is driven from a 
spiral ordered phase into a spin liquid phase by varying 
some parameter in the Hamiltonian, then the spin liquid 
phase has fluctuating spiral order and deconfined spin- 
1/2 spinons. In this scenario the RVB phase is protected 
by the gap to break a singlet bond and the transition to 
order occurs when the gap closes and an ordered moment 
can form from condensed spinons (or pairs of spinons). 
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FIG. 9: (color online) Phase diagram of a quasi-2D frustrated 
quantum magnet with deconfined spinons near an instability 
to spiral long-range order driven by a small parameter V in 
the Hamiltonian (such as the inter-layer coupling J"). The 
vertical dashed arrow shows the region that would be probed 
by experiments on CsgCuCU. Dashed lines mark cross-over 
to critical behavior from gapped fractional spin liquid for V < 
Vc and from renormalized classical on the ordered side {V > 
Vc)- At temperatures larger than the exchange energies, there 
is a cross-over to paramagnetic behavior (wavy dashed line). 



The magnetic order is an incommensurate spiral to min- 
imize the frustrated coupUngs and symmetry breaking in 
the spiral plane predicts a gapless Goldstone spin-wave 
mode at low energies, which can also be regarded as a 
two-spinon bound state stabilized by the weak mean-field 
effects of the magnetic order. In this phase diagram sce- 
nario deconfined spinons are expected to occur at higher 
energies above the scale of the binding potential, and at 
all energy scales in the spin liquid phase above Tn when 
the mean-field binding effects are cancelled by small ther- 
mal fluctuations, in agreement with experiments. 

Two theories of 2D spin liquid phases with deconfined 
spinons have so far been proposed for the main (J, J') 
Hamiltonian of CS2CUCI4, a model with gapped, bosonic 
spinonsii and a model with gapless, fermionic spinonsi^. 
In the latter case the phase diagram in Fig. |51 would be 
modified such that the fractional spin-liquid phase exists 
only at the critical point, i.e. Vc = 0. Another interest- 
ing theoretical possibility is that of a quantum phase with 
magnetic order and gapless magnons co- existing with 
topological order and gapped spinons recently discussed 
in the context of an unfrustrated square lattica^^. 



2. Paramagnetic scattering 

Neutron scattering studies of the high-temperature 
paramagnetic scattering are highly uncommon, whereas 



in fact they can reveal important information about the 
system. Here we consider a physical picture for the ob- 
served scattering in the paramagnetic phase by making 
an analogy with two model systems where a number of 
exact theoretical results are known. We refer in some 
detail to the ID Haldane-Shastry model (equivalent to 
an ideal ID spinon gas) and the XY chain (equivalent to 
a non-interacting lattice fermion gas) where the physics 
at T = cx) can be understood in terms of a disordered 
gas of quasiparticles where as a consequence of spin con- 
servation rules a neutron scattering process involves only 
two-particle states. 

The ID Heisenberg model with 1/r^ exchange 
(Haldane-Shastry model-HSM)i is equivalent to an ideal 
spinon gas. At T=0 the ground state is the spinon vac- 
uum and the spin-1 excitations (observable by neutron 
scattering) are two-spinon creation precessesi. At T = 00 
exact diagonalizations on finite clustergS^ find that the 
spin-1 excitations form a continuum in (fc,w) that ex- 
tends up to the same upper boundary as the two- 
spinon continuum at T=0; the scattering has a sharp 
cutoff at this boundary and no other states at higher 
energies are accessible. Combining this result with the 
fact that the Hilbert space of the HSM Hamiltonian can 
be entirely described in a spinon basis (all eigenstates 
can be written in terms of states with an even number 
of spinonsi) suggests that the dynamical correlations at 
T = 00 may also be physically interpreted in terms of 
two-spinon scattering processes. For example, S ^(k,uj) 
contains scattering events where (i) two up-spinons are 
created, (ii) two down-spinons are annihilated or (iii) an 
up-spinon is created and a down-spinon is annihilated; 
the intensity is a sum of the cross-sections for the above 
processes, weighted by the statistical probability of all 
possible initial states. Since in the HSM model the spinon 
dispersion is not affected by the presence of other spinons 
in the ground statei, the maximum two-spinon energy at 
T = 00 (when the ground state can be regarded as a 
dense gas of spinons) is the same as at T=0 (when the 
ground state is a spinon vacuum) , in agreement with the 
exact diagonalization results'^^ showing an upper bound- 
ary unaffected by temperature. 

Another system where only two-particle states con- 
tribute to dynamical correlations (of a conserved op- 
erator) is the XY chain H = JY.rSfSf+^ + S'fSf+i, 
which can be recast into a problem of non-interacting lat- 
tice fermions using a Jordan- Wigner spin-particle map- 
ping (up spin = occupied state, down spin — empty 
state). The Hilbert space consists of eigenstates with 
a fixed number of particles and the ground state at T=Q 
has all negative energy states occupied and all posi- 
tive energy states empty, where the single-particle dis- 
persion is ujk — Jcos(27rfc). The excitations contribut- 
ing to the dynamical correlations of the conserved 
operator consist of particle-hole pairgSiSi: S^^{k,uj) = 
Efci,fc2(l~/'=i)/'=2'^('^-'^fci+'^A;2)^(fc + fci-fc2)/A^, where 
fk = l/[exp{uJk/kBT) + 1] is the statistical probability 
that the state k is occupied. At T = 00 every state is 
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equally probable to be occupied or empty {fk = 1/2, the 
ground state is a disordered gas of fermions) and the S^^ 
dynamical correlations are a particle-hole pair continuum 
contained within — < w < with the same upper 
boundary as at T=0, ojj! = 2J\ sin(7rA:)|. 

The above interpretation of the paramagnetic scatter- 
ing in the ideal spinon gas (HSM model) in terms of two- 
spinon states could also be extended to the ID nearest- 
neighbor HAFC chain. In this case however interactions 
between spinous renormalize downwards the dispersion 
relation^. Therefore the upper boundary of the two- 
spinon continuum at T = oo may be slightly reduced 
compared to T=0 as a consequence of energy renormal- 
ization effects. This conjecture is consistent with exact 
diagonalizations on finite clusters^, which found that 
most of the scattering weight (within a few percent) at 
T = oo is located below the T—0 two-spinon boundary; 
the same study?? also found that at T = oo the lineshapes 
near the antiferromagnetic zone center approached a dis- 
tinctly non-Gaussian shape, better approximated by a 
plateau-like function centered at zero energy (similar in 
shape to the solid curve in Fig. |HJ but where the hor- 
izontal axis is energy). The lineshape observed in Fig. 
IHl^top right) is better described by a plateau-like func- 
tion (thick dashed line) than by a Gaussian of the same 
second moment (dash-dotted line). In this comparison 
both curves are centered at co—O and are multiplied by 
(cij/fcBr)/[l — exp(— CLi/fcsT)] to account for thermal pop- 
ulation effects to first order in Lu/kBT. 

The observation of continuum scattering extending up 
to an upper boundary that is only slightly reduced com- 
pared to the low-temperature extent of the continuum 
scattering [see Fig. IS^bottom)] suggests that the para- 
magnetic lineshapes in CS2CUCI4 may be understood in 
terms of (mainly) two-spinon scattering processes. 



VI. CONCLUSIONS 

In conclusion we have measured the magnetic excita- 
tions in the quasi-2D frustrated quantum antiferromag- 
net CS2CUCI4 using high-resolution time-of- flight neu- 
tron spectroscopy. In the 2D spin liquid phase above the 
magnetic ordering transition the dynamical correlations 
are dominated by highly-dispersive scattering continua, 
characteristic of fractionalization of spin waves into pairs 
of deconfined, S'=l/2 spinous. Boundaries of the scatter- 
ing continua are strongly-2D dispersive indicating that 
2D effects are important up to all energy scales of the 
excitations. Upon cooling into the 3D spiral ordered 
phase sharp magnon modes emerge at the lowest energies, 
but the dispersion relations and continuum scattering at 
medium to high energies were essentially unchanged com- 
pared to above Tjv. This suggests a dimensional crossover 
from 3D S=l spin waves at low energies to 2D S—1/2 
spinous at medium to high energies. A phase diagram 
scenario locating CS2CUCI4 in the proximity of a spinon 
de-confinement transition was proposed. 



We hope these results will stimulate further theo- 
retical work in the field of frustrated quantum mag- 
nets and fractionalized phases; knowing the Hamiltonian 
rigorously-"^^ should provide a test of the various theo- 
retical approaches. We also note that crucial in our ex- 
perimental exploration of the physics and in building up 
a phenomenological picture of the excitations was the 
use of time-of-flight neutron spectroscopy, which allowed 
probing the dynamical correlations over a very large part 
of phase space. 
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APPENDIX A: LINEAR SPIN- WAVE THEORY 
FOR A SPIRAL-ORDERED MAGNET 

In this section we describe the ground state and exci- 
tations of the frustrated Hamiltonian on the anisotropic 
triangular lattice in Fig.^b) using linear spin-wave the- 
ory. This approach assumes that spins order in a well- 
defined structure (found by minimizing the mean-field en- 
ergy) and treats excitations as local spin deviations that 
propagate by inter-site hopping. Linearizing the equa- 
tions of motion results in wave-like propagating modes, 
the magnon quasi-particles, which have quantized spin 
angular momentum S—1 and Bose statistics. We dis- 
cuss the dynamical correlations parameterized in terms 
of one- and two-magnon scattering processes. 

To find the mean-field ground state it is illuminating to 
consider first an isosceles triangle [see Fig. llUf a')]. which is 
the building block of the full two-dimensional anisotropic 
triangular lattice in Fig.^b). In the absence of the zig- 
zag coupling J'=0 the strong antiferromagnetic exchange 
J favors antiparallel alignment of the spins 2 and 3. This 
coUinear order becomes unstable in the presence of a cou- 
pling J' and the mean field created by spin 1 induces a 
spin-flop-like tilt of the other two spins as indicated in 
the figure. The angular tilt -0 depends on the relative 
strength J'/ J and is therefore incommensurate with the 
lattice. Generalization of this model to the full 2D lattice 
results in the non-coUincar spiral order indicated in Fig. 
llUf b). The ordered spin at site R is (Sr) — {S) cos(Q • 
R + <j)o)x + (S) sm{Q ■ R + (j)o)y, where (S) is the magni- 
tude of the ordered spin moment and 0o is an arbitrary 
phase angle, x and y are orthogonal unit vectors defin- 
ing the plane of spin rotation. The ordering wavevector 
Q is found by minimizing the exchange energy (per spin) 
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Jfc = Jcos(27rfc) + 2J'cos(7r/c) cos(7rZ), k = {h,k,l) and is 
Q = (0,0.5 + ec, 0) where = sm'^lJ' / {2J)]/tt is the in- 
commensuration relative to Neel order. In the mean-field 
approximation the incommensurate spiral order occurs 
immediately as the chains along J become coupled, i.e. 
J' 7^0, and is stable throughout the range < J' /J < 2, 
which includes the isotropic triangular lattice J' /J~l for 
which Ec = 1/6- 

An incommensurate spiral order as depicted in Fig. 
[TUT b) is observedi^ in CS2CUCI4 at temperatures be- 
low r;v=0.62(l) K. The ordered spins rotate in cy- 
cloids, nearly contained in the (6c) plane due to a small 
anisotropy. The angle between consecutive spins along 
& is 6* = 27rQ=180°-hl0.8(8)° which give an incommen- 
suration eo=0. 030(2). The measured incommensuration 
is significantly smaller than the classical value ec=0.054 
corresponding to the measured exchange parametersii 
J'/J=0.34(3). This large renormalization eo/ec=0.56 is 
due to quantum fluctuations not captured by the clas- 
sical mean-field approximation, but well reproduced by 
calculations which include nearest-neighbor spin-singlet 
correlations in the incommensurate ground statei^*^. 
Those correlations favor antiferromagnetic alignment of 
the paired spins, therefore acting to reduce the incom- 
mensuration as observed. 

Since the ordered spin direction rotates as one moves 
from site to site it is convenient to define a local ref- 
erence frame (f , 77, C) such that the mean spin direction 
at each site appears along the C a-xis for all sites. This 
transformation is illustrated in Fig. [rUT cl and is defined 
by 



COS ( 
sin (J 
5", 



- sine 
5*^ cos ( 



(Al) 



where the rotation angle at site i? is = Q • i?-|-0o. Spin 
deviation operators are expressed in terms of magnon 
Bose operators using the Holstein-Primakoff— formalism 



Sr = 



Si 
4 



i5^ = 



2Sa 



R 



^R 



S - ajfjflR, 



(A2) 



where . . . stand for higher-order terms involving three or 
more operators, creates a spin deviation at site R and 
its Fourier transform aj[ = J^r alie~^'''^/\/N creates an 
extended wave {N is the total number of spins). 

In a 1/ S expansion the leading term in the Hamiltonian 
becomes a quadratic form of magnon operators. This can 
be diagonalized using standard techniques to obtain the 
wave- vector-dependent spin-wave energies^ 



Uk = 2S^{Jk - Jq) [{Jk-Q + Jk+q) 12 - Jq], (A3) 

where Jfc — Jcos(27rfc) -I- 2 J' cos(7rfc) cos(7r/) is the 
Fourier transform of the exchange couplings. Higher- 
order terms in the Hamiltonian amount to interactions 
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FIG. 10; (a) Schematic diagram of the mean-field ground 
state for three spins located in the corners of an isosceles 
triangle with antiferromagnetic couplings J (horizontal thick 
bond) and J' (thin zig-zag bonds), (b) Spiral order on the 
anisotropic triangular lattice obtained by generalizing the pic- 
ture shown in (a). Thick arrows indicate the ordered spin 
direction, which rotates in the {x,y) plane by an incommen- 
surate angle 6 — 2nQ between neighboring sites in the direc- 
tion of the strong exchange J (horizontal thick bonds), (c) 
Definition of the rotating reference frame {^,r],Q where ( is 
along the local mean spin direction and rj is perpendicular to 
the rotation plane. 



between magnons, which are expected to produce (rel- 
atively) small renormalizations of the above spin-wave 
energies. A typical plot of the dispersion relation in Eq. 
(|A3p is shown in Fig. Ela). The dispersion has a sinu- 
soidal shape with gapless modes at the ferromagnetic 
zone center fe = and at the incommensurate order- 
ing points k — zLQ in the 2D Brillouin zone. The luq 
mode corresponds to a uniform rigid rotation of all spins 
in the ordering plane and is called a phason, whereas the 
^±Q modes correspond to fluctuations where the plane 
in which spins rotate cants away from the ordering plane. 
Those three gapless excitations are Goldstone modes as- 
sociated with the breaking of symmetry in spin space by 
the spiral long-range order [in the presence of a small 
easy-plane anisotropy the mode uj±q acquires a gap but 
uq remains gapless as it involves spin rotations inside the 
easy-plane] . 

Information about the nature and properties of the 
spin excitations is contained in the dynamical correla- 
tion function 3°^^ {k,Lu), directly measured by neutron 
scatteringSS,. This is defined as the Fourier transform of 
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the space and time correlation function 

S-f'ik^u;) = — / dtJ2{SS{0)S'^{t))e 

(A4) 

where a and (3 label Cartesian axes and (. . .) denotes 
a ground state average. Because the Hamiltonian 
conserves the total spin component S^, the ofF-diagonal 
terms cancel out S°'^{k, w) = for a 7^ /3, see Ref. l27l 

One-magnon (lA^) excitations occur in the fluctua- 
tions polarized transverse to the ordered spin direction 

The transverse dynamical correlations at T=0 are ob- 
tained as^l 



S 



Sl%^{k,uJ)^— \uk + Vk\^ 5{uj - ujk), 
S 



(A5) 



Sharp one-magnon peaks occur at wavevectors and ener- 
gies given by the dispersion relation ojk. The scattering 
intensities depend on the functions Uk and Vk defined as 
Uk — cosh 9k and Vk = sinh 6^ where tanh 29^ = Bf, / Ak 
and 

Ak = 25 + - [Jfc_Q + Jfc+Q] - Jq| , 

Bk = 25 — - [Jje-Q + Jfe+g] 



The dispersion relation in Eq 
expressed as uJk 



SJ can alternatively be 



Fluctuations in the longitudinal spin component (|| C) 
are expressed as a product of two local magnon oper- 
ators in Eq. l|A2p . After transformation to normal op- 
erators that diagonalize the quadratic Hamiltonian, the 
longitudinal dynamical correlation function is expressed 
in terms of two-magnon (2A^) scattering processes. At 
r=0 only two-magnon creation processes can occur and 
their cross-section is given 



xS{uo - uoki ~ WfcJJ(fc + ki 



k. 



(A6) 



where /(fci, ^2) = \u-k1Vk2 + Uk^V-ki \ is the scattering 
weight. This formula is valid in the harmonic approxima- 
tion where magnon-magnon interactions are neglected, 
i.e. only the quadratic part is retained in the Hamiltonian 
expansion in terms of magnon operators. In Eq. ljA6p ki 
and k2 are the wave vectors of the two magnons, which 
can be created anywhere in the first Brillouin zone and r 
is a vector of the reciprocal lattice. For a given wavevec- 
tor transfer k two-magnon processes contribute an ex- 
tended scattering continuum at energies above a lower 
threshold. This threshold corresponds to creating one 
magnon in the lowest energy state (in this case u;fci=0 at 
a wave vector ki = 0, —Q or +Q) and the other magnon 
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FIG. 11: (a) Two-magnon scattering intensity ^^^(fc — <5,a;) 
in Eq. JXgJ as a function of energy and wave vector along b* . 
Density of scattered points represents intensity. The plotted 
points are random two-magnon scattering events generated 
using a Monte Carlo algorithm for the sum in Eq. lElJ 

as 

described in the text. The displayed two-magnon scatter- 
ing continuum is bounded at low energies by the lower of 
the three dispersion relations LOk (solid line), L0k-2Q (dashed 
line) and cuk-Q (dash-dotted line), (b) Energy-integrated 
two-magnon scattering intensity as a function of wave vec- 
tor along b* . Solid points represent J dujS2j^{k — Q,u)) 
calculated using a distribution of scattering events such as 
shown in (a). For comparison, the intensity of the in-plane 
magnon with dispersion Uk-q is shown by the dash-dotted 
line (/ dujSl%f{k-Q ,uj) see Eq. lASH . Calculations are plot- 
ted for renormalized couplings J=0.61 meV and J' / J— 0.175 
that fit the observed dispersion relation in CS2CUCI4 plotted 
in Fig. 13 One- and two-magnon scattering intensities were 
normalized against sum rules as described in the text with 
AS/S ~0.25 estimated from experimental-; the term {uv}^ 
was calculated to be very small and was neglected. 



at another place in the zone. The sum in (jA6p was calcu- 
lated numerically using a weighted Monte Carlo method 
described in Ref. [s^ From a large number of randomly- 
generated two-magnon scattering processes a distribution 
of 2 X 10^ events was selected according to a probabil- 
ity proportional to the scattering weight /(fei,fe2); this 
gave the complete picture of the two-magnon intensity 
in the whole Brillouin zone. Fig. Illf a) shows a plot of 
the resulting intensity shifted in wave vector by Q to ob- 
tain the quantity 52^(fe — Q, oj) relevant for comparison 



16 



with experiments. The intensities obtained after numeric 
summation were then normaUzed to satisfy the sum rule 
for the total two-magnon scattering in the Brillouin zone, 
as discussed below. 

Longitudinal inelastic processes described in Eq. (|A6p 
can occur because the spin moment is not fully ordered 
in the ground state. This disordering is due to zero- 
point quantum fluctuations and the spin reduction cal- 
culated in the linear spin-wave approximation is AS* = 
Sfe bfcP/-^ where k spans the Brillouin zone. Applied 
to the relevant 2D coupHngs in CS2CUCI4 ( J'/ J=0.34(3)) 
this gives a very large spin reduction A5'=0.43, see 
Ref,40,, however mean- field effects from the inter-layer 
couplings and anisotropy terms not included in this cal- 
culation can quench some of the low-energy fluctua- 
tions and partially restore order, leading to a smaller 
AS*. The total two-magnon scattering intensity inte- 
grated over energy and wavevector in a Brillouin zone is 
AS{l + AS) + {uvy (per spin), where (uv) = J2k UkVk/N 
is in general non zero if the ordering wavevector Q is 
incommensurate. Zero-point fluctuations also affect the 
one-magnon scattering by reducing the intensity factor 



in Eq. from S to S = S - AS ^ 



[uv) 



\l + 2ASy 



With this normalization the sum rule for the total scat- 
tering in the Brillouin zone S{S + 1) is exhausted by 
elastic Bragg scattering [S — ASY shared between the 
two positions k = ±Q, inelastic one-magnon 5'(1-|-2A5') 
and two-magnon processes. Higher-order scattering pro- 
cesses involving three or more magnons are neglected here 
since their scattering weight decreases very rapidly with 
increasing particle number. Those processes redistribute 
some of the scattering weight in the zone towards higher 
energies. 

To obtain the dynamical correlations measured in an 
experiment, Eqs. ljA5|l and ljA6|l are converted to the 
fixed laboratory reference frame (x, y, z) using the trans- 
formation in Eq. (|A1|) . Sharp one-magnon peaks occur 
both in the out-of-plane (|| z) 



5{u~uk), (A7) 



2 uj 

as well as in the in-plane (|| x,y) dynamical correlations 



_ 1 
4 



Sli,{k-Q,uj) + Sli^{k + Q,Lo 



S Ak-Q + Bk-Q 



S Ak+Q 
8 



k+Q 



- LUk-q) + 



5{uj-ujk+Q) (A8) 



In total three spin-wave modes are observed: a principal 
mode ujk polarized out-of-plane {\\ z) and two secondary 
modes cj^ — ujk-Q and cj^ = ^k+Q polarized in-plane 
(II x,y). Those three modes can also be regarded as ex- 
citations with quantized spin moment along the z-axis 




FIG. 12: (color online) Contour plot of the 2D dispersion rela- 
tion using a square-lattice representation of Fig.^b) (J' along 
the square edges). Filled circles are zone centers and stars 
show positions of incommensurate magnetic Bragg peaks. 



normal to the spiral plane, i.e. the principal mode LUk has 
spin S^=0 and the two secondary modes uj^ have spin 
5^ = ±1 {+z defines the sense of rotation in the spiral). 
The three dispersions are plotted along the b* and c* di- 
rections in Fig.Oand for all modes the intensity decreases 
as 1 /uj with increasing energy. The secondary modes are 
images of the main mode displaced in wavevector by — Q 
and +Q. This apparent splitting into three modes is a 
characteristic feature of systems with helical order, i.e. 
helix, cycloid, and cone^. 

The two-magnon scattering is polarized in-plane and 
is obtained using transformation (|A1|) as 



S^X,ik,io) 



SiUk-Q.oj) 



(A9) 



Starting from the two-magnon continuum S'2^(fc,w) in 
Eq. ljA6|l calculated in the rotating reference frame the 
wave vector values k are shifted by -fQ and —Q- This 
gives two continua, their superposition is the two-magnon 
scattering in the fixed reference frame in Eq. l|A9p . Fig. 
Illf a) shows a plot of the first term in Eq. ljA9|l as a 
function of energy w and wave-vector k along the b* di- 
rection. In the figure the density of points indicates scat- 
tering intensity, which decreases rapidly with increasing 
energy. The energy-integrated two-magnon scattering in- 
tensity is plotted in Fig. Ilir bVsolid points): the inten- 
sity is largest near the magnetic Bragg peak position and 
cancels at the ferromagnetic zone center fc — > r. For 
5=1/2 and AS'=0.125 the two-magnon spectral weight is 
throughout most of the Brillouin zone comparable to the 
weight of the in-plane polarized one-magnon mode to^ 
(dashed-dotted line) that occurs near the lower boundary 
of the two-magnon continuum, whereas for wavevectors 
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near the magnetic Bragg peak position most scattering 
weight is in the one-magnon channel. 

Including the polarization factors for magnetic scatter- 
ing the one-magnon cross-section measured by neutron 
scattering is 

hM{k,uj)=p^S^,Xiik,^)+P.SrMik,^) (AlO) 
and the two-magnon cross-section is 

l2M{k,uj) ^p^S^^{k,uj) (All) 
where the polarization factors are 

Px ^ ^ + cos^ ak 

p^ = sin^ttfc (A12) 

The relative intensity of in-plane (|| x) and out-of-plane 
(II z) polarized excitations can therefore be changed by 



varying the angle ak between the scattering wave vec- 
tor k and the axis z normal to the rotation plane. For 
CS2CUCI4 the ordered spins rotate nearly in the (be) 
plane, so the axes are identified as {x,y,z) = (6, c, a). 
For comparison with experiments the magnetic form fac- 
tor |/(A;)p for Cu^+ ions^^ was also included on the r.h.s. 
in Eqs. (|IlO|l and ifnTjl . 



APPENDIX A: SQUARE-LATTICE NOTATION 

Theoretical models^**'^^ for the triangular lattice in Fig. 
n^b) often regard the problem as a square-lattice with 
exchange J' is along the edges and J along one of the 
diagonals. A contour plot of the 2D dispersion relation 
in Eq. ^ using the square-lattice basiSi^S is shown in Fig. 

[El 



^ F.D.M. Haldane and M.R. Zirnbauer, Phys. Rev. Lett. 71, 
4055 (1993). 

2 G. Miiller, H. Thomas, H. Beck and J.C. Bonner, Phys. 
Rev. B 24, 1429 (1981); 

3 D.A. Tennant, R.A. Cowley, S.E. Nagler and A.M. Tsvelik, 
Phys. Rev. B 52, 13368 (1995); D.C. Dender, D. Davidovic, 
D.H. Reich, C. Broholm, K. Lefmann, and G. Aeppli, Phys. 
Rev. B 53, 2583 (1996). 

* F.D.M. Haldane, J. Appl. Phys. 57, 3359 (1985). 
^ P.W. Anderson, Mat. Res. Bull. 8, 153 (1973). 

® S.M. Hayden, G. Aepph, R. Osborn, A.D. Taylor, T.G. 
Perring, S.W. Cheong, and Z. Fisk, Phys. Rev. Lett. 67, 
3622 (1991); H.M. R0nnow, D.F. McMorrow, R. Coldea, 
A. Harrison, I.D. Youngson, T.G. Perring, G. Aeppli, O. 
Syljuasen, K. Lefmann, and C. Rischel, Phys. Rev. Lett. 
87, 037202 (2001). 

H. Kageyama, K. Yoshimura, R. Stern, N. V. Mushnikov, 
K. Onizuka, M. Kato, K. Kosuge, C.P. Slichter, T. Goto, 
and Y. Ueda, Phys. Rev. Lett. 82, 3168 (1999). 

* R.H. McKenzie, Comment. Condens. Matter Phys. 18, 309 
(1998). 

® K. Takada, H. Sakurai, E. Takayama-Muromachi, F. 
Izumi, R.A. Dilanian, and T. Sasaki, Nature 422, 53 
(2003). 

° R. Coldea, D.A. Tennant, A.M. Tsvelik, and Z. Tylczynski, 

Phys. Rev. Lett. 86, 1335 (2001). 
^ V. Kalmeyer and R.B. Laughlin, Phys. Rev. Lett. 59, 

2095 (1987); S. Sachdev, Phys. Rev. B 45, 12377 (1992); 

R. Moessner and S.L. Sondhi, Phys. Rev. Lett. 86, 1881 

(2001); S. Sachdev and K. Park, Ann. Phys. New York 

298, 58 (2002); 
^ N. Read and S. Sachdev, Phys. Rev. Lett. 66, 1773 (1991). 
^ A.A. Nersesyan and A.M. Tsvelik,Phys. Rev. B 67, 024422 

(2003). 

" C.H. Chung, J.B. Marston and R.H. McKenzie, J. Phys.: 
Condens. Matter 13, 5159 (2001). C.H. Chung, K. Voelker 
and Y.B. Kim, cond-mat/0302452 (2003). 
^ Yi Zhou and X.G. Wen, cond-mat/0210662 
^ T. Senthil and M.P.A. Fisher, Phys. Rev. B 62, 7850 



(2000). 

R. Coldea, D.A. Tennant, K. Habicht, P. SmeibidI, C. 
Wolters and Z. Tylczynski, Phys. Rev. Lett. 88, 137203 
(2002). 

S. Bailleul, D. Svornos, P. Porcher and A. Tomas, C. R. 
Acad. Sci. Ser. II 313, 1149 (1991). 

R. Coldea, D.A. Tennant, R.A. Cowley, D.F. McMorrow, 

B. Dorner, and Z. Tylczynski, J. Phys.:Condens. Matter 8, 
7473 (1996). 

L.V. Soboleva, L.M. Belaev, V.V. Ogadzhanova and M.G. 
Vasileva, Krystallografiya 26, 817 (1981). 

C. J. Garble and M.A. Adams, Physica B 182, 431 (1992); 
R. Crevecoeur, I. de Schepper, L. de Graaf, W. Montfrooij, 
E. Svensson, and C. Carlile, Nucl. Instrum. Meth. A 356, 
415 (1995). 

The transformation of wavevector components from rlu 
units k = {h, k, I) to the square-lattice representation of 
Fig.^b) (J' along the square edges) is kx = {k — 1)-k and 
ky = {k-\-l)Ti (for unit square-lattice spacing). The Fourier 
transfer of the exchange couplings is = J cos(A;3: + ky)-\- 
J'{coskx + cos ky), the incommensurate ordering wavevec- 
tor is Q = (-1 -I- eoTT, f -I- eoTi") and "saddle" points of the 
2D dispersion relation in Eq. ||5J occur near (Jzx^ky) — 
(t"' T") + ("^'0"", with / an even integer. 
Experiments made with the sample aligned in the {a* ,b*) 
plane used a spectrometer configuration with slightly nar- 
rower resolution Ai?=0. 016(1) meV. 
R.R.P. Singh, Phys. Rev. B 39, 9760 (1989). 
International Tables for Crystallography, edited by A. J.C. 
Wilson (Kluwer, Dordrecht, 1995), Vol.C, p. 391. 

A. V. Chubukov, T. Senthil, and S. Sachdev, Phys. Rev. 
Lett. 72, 2089 (1994); A.V. Chubukov, S. Sachdev, and T. 
Senthil, Nuclear Phys. B 426, 601 (1994). 

S.W. Lovesey, in Theory of Neutron Scattering from Con- 
densed Matter (Clarendon, Oxford, England, 1987). 
P.G. de Gennes, J. Phys. Chem. Solids 4, 223 (1958). 

B. A. Bernevig, D. Giuliano, R.B. Laughlin, Phys. Rev. B 
64, 024425 (2001); 

M. Bocquet, F.H.L. Essler, A.M. Tsvelik, and A.O. 



18 



Gogolin, Phys. Rev. B 64, 094425 (2001); M. Bocquet, 
Phys. Rev. B 65, 184415 (2002). 

F.H.L. Essler, A.M. Tsvelik, and G. Delfino, Phys. Rev. B 

56, 11001 (1997). 

C. Lannert and M.P.A. Fisher, lcond-m at/0204382| 

K. Fabricius, U. Low, and J. Stolze, Phys. Rev. B 55, 5833 

(1997). 

S. Katsura, T. Horiguchi, and M. Suzuki, Physica 46, 67 
(1970). 

Zheng Weihong, R.H. McKenzie, and R.P. Singh, Phys. 
Rev. B 59, 14367 (1999). 

T. Holstein and H. Primakoff, Phys. Rev. 58, 1098 (1940). 
T. Nagamiya, Solid State Physics, edited by F. Seitz, D. 



Turnbull and H. Ehrenreich, (Academic, New York, 1967), 
Vol. 20, p. 305; J. Jensen and A.R. Mackintosh, Rare Earth 
Magnetism. Structures and Excitations, (Clarendon, Ox- 
ford, 1991), p. 286. 

I.U. Heilmann, J.K. Kjems, Y. Endoh, G.F. Reiter, G. 
Shhane, and R.J. Birgeneau, Phys. Rev. B 24, 3939 (1981). 
D.A. Tennant, S.E. Nagler, D. Welz, G. Shirane, and K. 
Yamada, Phys. Rev. B 52, 13381 (1995). 
J. Merino, R.H. McKenzie, J.B. Marston and C.H. Chung, 
J. Phys.: Condens. Matter 11, 2965 (1999); A.E. Trumper, 
Phys. Rev. B 60, 2987 (1999). 



